Differential-geometric approach to the integrability of 
hydrodynamic chains: the Haantjes tensor 



E.V. Ferapontov and D.G. Marshall 



Department of Mathematical Sciences 
Loughborough University 
Loughborough, Leicestershire LE11 3TU 
United Kingdom 

e-mails: 

E.V. FerapontovOlboro .ac.uk 
D.G. Marshall@lboro .ac.uk 



Abstract 

The integrability of an m-component system of hydrodynamic type, u t = V(u)u x , by 
the generalized hodograph method requires the diagonalizability of the m x m matrix U(u). 
This condition is known to be equivalent to the vanishing of the corresponding Haantjes 
tensor. We generalize this approach to hydrodynamic chains — infinite-component systems of 
hydrodynamic type for which the ooxoo matrix V(u) is 'sufficiently sparse'. For such systems 
the Haantjes tensor is well-defined, and the calculation of its components involves finite 
summations only. We illustrate our approach by classifying broad classes of conservative 
and Hamiltonian hydrodynamic chains with the zero Haantjes tensor. We prove that the 
vanishing of the Haantjes tensor is a necessary condition for a hydrodynamic chain to possess 
an infinity of semi-Hamiltonian hydrodynamic reductions, thus providing an easy-to-verify 
necessary condition for the integrability. 

MSC: 35L40, 35L65, 37K10. 
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1 Introduction 



Hydrodynamic chains are quasilinear first order PDEs of the form 

u t = V(u)u x (1) 

where u = (u 1 ,^ 2 , ...)* is an infinite-component column vector and V(u) is an oo x oo matrix. 
A classical example is the Benney chain (or Benney's moment equations), 

< = <+ 1 + (n-lK- 1 ^, (2) 

n = 1,2,..., which was derived in [2] from the equations for long nonlinear waves on a free 
surface. It was observed later in [52} ITT] that the same system results from a kinetic Vlasov 
equation. The system ([2]) has been thoroughly investigated in the subsequent publications 
[55], 123 (28J [331 E21 ESl 09J [T7] where, in particular, its Hamiltonian and integrability aspects 
were uncovered. Hydrodynamic reductions of the chain ([2]) were studied in |T8l [EH EUJ [26], [36] . 
Various deformations of Benney's equations are known. These include the modified Benney 
chain, 

u n t = u n x +1 + u l u n x + (n - l)u n ul, (3) 

obtained in [3UJ as a quasiclassical limit of the modified KP hierarchy. Its two-parameter defor- 
mation, 

n™ = u n x +1 + u x u n x + (a(n - 1) + b)u n u\, (4) 

was constructed in [29] along with further examples of Hamiltonian integrable chains possessing 
complete systems of commuting integrals. Another deformation scheme, based on the i?-matrix 
approach, was proposed in [3], see also [HIES]- The specialization of the chain corresponding 
to a = 0, b = —1, 

u n t = < +1 + u x u n x - u n u x x , (5) 

naturally appears in the theory of finite-gap solutions of integrable hierarchies of the KdV 
type [HE?]. Reductions of the chain both hydrodynamic and differential, were extensively 
investigated in [391 Wf[ 148] . The case a = 1, 6=2 arises from the kinetic model for rarefied 
bubbly flows [23] . see also [SU] for an alternative representation of this chain. 

A broad class of new examples was found in [3DJ, see also [1], based on the symmetry 
approach. These papers provide a classification of conservative chains of the form 

u}=ul, u\ = g{u x ,u 2 ,u s ) x , u\ = h{u x ,u 2 ,u s ,u A ') x ,..., (6) 

which are embedded into a commutative hierarchy of Egorov's type. It was observed that the 
function g(u l ,u 2 ,n 3 ) uniquely determines all other equations of the chain ([6]), as well as the 
whole associated hierarchy. Moreover, the function giu 1 ,u 2 ,u 3 ) was shown to satisfy an over- 
determined involutive system of third order PDEs whose general solution was expressed in terms 
of theta functions and solutions to the Chazy equation. 

Our approach to the integrability of hydrodynamic chains is motivated by the theory of 
finite-component systems of hydrodynamic type, that is, equations of the form ([1]) where u = 
(u , u 2 , u m )* is an m-component column vector and V^(u) = Vj(u) is an m x m matrix. 
Explicitly, one has 

ut = v j( u ) u li i,j = l,...,m; (7) 
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here and below we adopt the standard summation convention over the repeated upper and lower 
indices. Such systems naturally occur in applications in gas dynamics, fluid mechanics, chemical 
kinetics, Whitham averaging procedure, differential geometry and topological field theory. We 
refer to [511 l6l IT] for a further discussion and references. In what follows we consider the strictly 
hyperbolic case when the eigenvalues of the matrix Vj, called the characteristic speeds of the 
system (|?_J), are real and distinct. 

It is known that many particularly important systems of the form ([7]) are diagonalizable, 
that is, reducible to the Riemann invariant form 

Rl = X\K)Rl (8) 

where the characteristic speeds A 2 (R) satisfy the so-called semi-Hamiltonian property [51], 



dk = d/dR k , i ^ j ^ k. We emphasize that the semi-Hamiltonian property ([9]) is usually auto- 
matically satisfied for diagonalizable systems of the 'physical' origin. For instance, a conservative 
diagonalizable system is necessarily semi-Hamiltonian (see Appendix 1 for more details). Such 
systems possess infinitely many conservation laws and commuting flows of hydrodynamic type 
and can be linearized by the generalized hodograph method |51| . Their analytic, differential- 
geometric and Hamiltonian aspects are well-understood by now. 

Remarkably, there exists an efficient tensor criterion of the diagonalizability which does not 
require the actual computation of eigenvalues and eigenvectors of the matrix Vj. Let us first 
calculate the Nijenhuis tensor of the matrix Vj, 

N] k = v?d uP vl - v\d uV v) - v^vl - d^v?), (10) 
and introduce the Haantjes tensor 

H] k = N;y jV i - Nj r VpVf. - iv>x + N » r P - (ii) 

For strictly hyperbolic systems the condition of diagonalizability is given by the following theo- 
rem. 

Theorem 1 [22] A hydrodynamic type system with mutually distinct characteristic speeds is 
diagonalizable if and only if the corresponding Haantjes tensor ( T77]) is identically zero. 

Since components of the Haantjes tensor can be calculated (using computer algebra) in any 
coordinate system, this provides an efficient diagonalizability criterion. 



Our main observation is that both tensors (I10p and (llip make perfect sense for infinite 
matrices which are 'sufficiently sparse'. To be more rigorous, let us give the following 

Definition 1. An infinite matrix V(u) is said to belong to the class C (chain class) if it satisfies 
the following two natural properties: 

(a) each row ofV(u) contains finitely many nonzero elements; 

(b) each matrix element ofV(u) depends on finitely many variables u % . 

Notice that the chains (EJ) - ((6j) clearly belong to the class C. For matrices from the class C all 
contractions in the expressions (|10p and (jlip reduce to finite summations so that each particular 
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component H l - k is a well-defined object which can be effectively computed. Moreover, for a fixed 
value of an upper index i there exist only finitely many non-zero components H l - k . 
We propose the following 

Definition 2. A hydrodynamic chain from the class C is said to be diagonalizable if all compo- 
nents of the corresponding Haantjes tensor ill]) are zero. 

We point out that the chains © - © ar e diagonalizable in this sense. As we prove in Sect. 5, 
the vanishing of the Haantjes tensor is a necessary (and in many particularly important cases — 
sufficient) condition for a hydrodynamic chain to possess an infinity of finite-component diag- 
onalizable hydrodynamic reductions, thus justifying the above definition. The main advantage 
of our approach is its 'intrinsic' character: it does not require the knowledge of any 'extrinsic' 
objects such as commuting flows, Hamiltonian structures, Lax pairs, etc. 

The vanishing of the Haantjes tensor turns out to be an efficient classification criterion. As 
an elementary example let us consider the chain 

v,™ = + v^u™ + c n u n u\ 

where c n = const. One can readily verify that the vanishing of the Haantjes tensor implies the 
recurrence relation c n+ 2 = 2c n+ i — c n . Setting ci = b, C2 = a + b we recover the integrable chain 

Based on the same criterion, in Sect. 2 we classify diagonalizable chains of the type ([6]). 
It turns out that the conditions Hj k = are already sufficiently restrictive and imply an over- 
determined system expressing all second order partial derivatives of the function h in terms of g, 
see (I14D . The consistency conditions of these equations lead to a closed- form involutive system 
expressing all third order partial derivatives of g in terms of its lower order derivatives, see (|15j) 
(we have extensively used computer algebra to calculate the Haantjes tensor and to verify the 
involutivity by calculating the compatibility conditions). The requirement of the vanishing of 
other components H* k , i > 2, imposes no additional constraints on h and g: these conditions 
reconstruct the remaining equations of the chain ([6]). For instance, the conditions H 2 k = 
specify the right hand side of the fourth equation u\ = etc. 

The characterization of diagonalizable chains of a more general form, 

u\ = f(u x ,u 2 ) x , u\ = g(u l ,u 2 ,u 3 ) x , u\ = h(u l , u 2 , u 3 , u 4 ) x , (12) 

etc, is proposed in Sect. 3. As in the previous example, the conditions Hj k = lead to 
expressions for all second order partial derivatives of h in terms of g and /. The consistency 
conditions of these equations result in an involutive system expressing all third order partial 
derivatives of g and / in terms of lower order derivatives thereof. 

In Sect. 4 we classify diagonalizable Hamiltonian chains of the form 

* \ dx dx ) <9u' ^ ^ 

here B l i = (i — l)uj_|_ ? _2 and ^(u 1 , u 2 , -u 3 ) is a Hamiltonian density. The Benney chain 
corresponds to h = (n 3 + (u 1 ) 2 )/2, see |28| . We have found a broad family of new integrable 
Hamiltonian densities, thus extending the results of |29j. 

In Sect. 5 we prove that the condition of diagonalizability is necessary for the existence 
of 'sufficiently many' hydrodynamic reductions. Recall that an m-component hydrodynamic 
reduction of an infinite chain is specified by parametric equations 



u 



1 = u 1 (R\...,R m ), u 2 = u 2 {R\...,R m ), u 3 = u 3 (R\...,R ri 
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etc, where the Riemann invariants R 1 ,... ,R m solve a diagonal system (|H|) whose characteristic 
speeds satisfy the semi-Hamiltonian property 0. It is required that all equations of the chain 
are satisfied identically modulo ([8]). Thus, an infinite chain reduces to a system with finitely 
many dependent variables. It was demonstrated in [\9\ that the Benney chain ([2]) possesses 
infinitely many m-component reductions of this type parametrized by m arbitrary functions of 
a single variable. The same is true for other chains ([2]) - ©. Based on these and other examples 
we propose the following 

Definition 3. A hydrodynamic chain from the class C is said to be integrable if, for any 
m, it possesses infinitely many m-component semi-Hamiltonian reductions parametrised by m 
arbitrary functions of a single variable. 

Remark. It was observed in [32^ 02] that Hamiltonian chains (I13p possess m-component re- 
ductions for any Hamiltonian density h, even in the non-integrable case. The crucial difference 
is that for integrable chains these reductions are semi-Hamiltonian and depend on m arbitrary 
functions of a single variable, while in the non-integrable situation for any m there exists a 
unique m-component reduction which is not diagonalizable. 

In Sect. 5 we prove our main result: 

Theorem 2 The vanishing of the Haantjes tensor H is a necessary condition for the integra- 
bility of hydrodynamic chains from the class C . 

If the spectrum of the infinite matrix V is simple, that is, for a generic A there exists a unique 
formal eigenvector £ such that V£, = A£ (we point out that all of the above examples satisfy this 
property), one has the following stronger result: 

Theorem 3 In the simple spectrum case the vanishing of the Haantjes tensor H is necessary and 
sufficient for the existence of two-component reductions parametrized by two arbitrary functions 
of a single variable. 

Theorem 2 provides an easy-to-verify necessary condition for testing the integrability of 
hydrodynamic chains. We emphasize that the vanishing of the Haantjes tensor is not sufficient 
for the integrability in general: one can construct examples of diagonalizable chains which 
possess infinitely many diagonal reductions none of which are semi-Hamiltonian (see Sect. 5). 
To eliminate these pathological cases let us recall that for finite-component systems ([7]) there 
exists a tensor object which is responsible for the semi-Hamiltonian property [H]. This is a 
(l,3)-tensor P^- (see Appendix 1 for explicit formulas in terms of the matrix v 1 -). Similarly to 
the Haantjes tensor H, the tensor P is well-defined for hydrodynamic chains from the class C. 
We conclude this introduction by formulating the following 

Conjecture. The vanishing of both tensors H and P is necessary and sufficient for the inte- 
grability of hydrodynamic chains from the class C . 

The necessity of this conjecture (that is, the statement that the integrability implies the vanishing 
of both H and P) is a relatively simple fact, see Sect. 5 for a proof. The sufficiency is a far 
more delicate property which we were not able to establish in general. We point out that the 
vanishing of H alone (in fact, the vanishing of the very few first components of H), is already 
sufficiently restrictive and implies the integrability in many particularly important cases (e.g. 
for conservative chains, Hamiltonian chains, etc). Recall that, for finite-component systems of 
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hydrodynamic type, our conjecture is a well-known result: any diagonalizable semi-Hamiltonian 
system possesses infinitely many conservation laws and commuting flows of hydrodynamic type, 
and can be solved by the generalized hodograph method [51J. 



2 Classification of diagonalizable chains of the type ([6]) 

The structure of equations ([6]) implies that the only nonzero components H l - k of the Haantjes 
tensor are the ones with j < i + 4, k < i + 4. Taking into account the skew-symmetry of the 
Haantjes tensor in its lower indices, this leaves ten essentially different components of the type 
Hj k . Equating them to zero we obtain the expressions for all of the ten second order partial 
derivatives of h(u , u 2 , u 3 , u 4 ): 



/in 
h 2 2 



2giffi2 ~ 92911 + 2/tiffi 3 

53 

92291 + 9U + 913^2 + 923^-1 
93 

92922 + 2gi2 + 2g 2 3^2 
93 

913(^3 - 92) + 92391 + 91293 + 933^1 



, 9l3\ n 3 — 92) -r 92391 -r 91293 -r 933"1 M 

"13 = > l i4 J 

93 

, .913 + h 3 g 2 3 + h 2 g33 

"23 — 922 H 



h-33 — 2^23 



93 

933(92 - 2/t 3 ) 

93 



, ^-4913 , h 4 g 2 3 , ft-4933 , „ 

"14 — , "24 — j "34 — , "44 — U. 

93 93 93 

Notice that these equations can be compactly written as 

2 1 

d 2 h = — (dhdg 3 + dgidu 2 + dg 2 dg - -g 2 d 2 g); 
93 2 

here both sides of the equality are understood as symmetric two-forms, and dg, dh, d 2 g, d 2 h 
denote the first and second symmetric differentials of g and h. The consistency conditions for 
the equations (fill) lead to closed-form expressions for all third order partial derivatives of the 
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function g(u ,u 2 ,u 3 ) in terms of its lower order derivatives: 

2^3 2513533 2523533 

5333 — j 5133 — j 5233 — > 

53 53 53 

25l3 2^13^23 

5113 — 5 5123 — , 5223 — > 

53 53 53 

2 

5222 = — (52523 + 523(53522 + 2^13) - 533(52522 + 2^12)) , 
53 

2 

5122 = — (51523 + 513(53522 + 513) - 533(51522 + 5n)) , (15) 
53 
2 

5112 = — (533(52511 - 251512) - 513(52513 - 253512) - 523(53511 - 251513)) , 
53 

5111 = — ((51 + 5l)5l3 + 5l 523 + 53(5l2 - 511522) - 52253351 

53 

+51353(511 + 2(51522 - 52512)) + 2523(52(53511 - 51513) - 5153512) 

- 533 ((51 + 5l)5ll - 25152512)) • 

This system is in involution and its general solution depends on ten integration constants, indeed, 
the values of 5 and its partial derivatives up to the second order can be prescribed arbitrarily at 
any fixed point Uq,Uq,Uq. The system (|15p was first derived in [ID] from the requirement that 
the chain © is embedded into a hierarchy of commuting hydrodynamic chains of Egorov's type. 
Exactly the same equations for 5 were obtained in |14j by applying the method of hydrodynamic 
reductions to the (2+l)-dimensional PDE 

Utt = g(Uxx,Uxt,Uxy) (16) 



which is naturally associated with the chain (J6J) ; here the function g is the same as in ([6]), (|15p . 
Thus, for hydrodynamic chains of the type ([6]) the condition of diagonalizability is necessary 
and sufficient for the integrability. 

One can show that the vanishing of other components of the Haantjes tensor does not 
impose any additional constraints on the derivatives of 5 and h. Thus, writing the fourth 
equation of the chain (J6j) in the form u\ = s(u 1 ,u 2 ,u 3 ,u' i ,u 5 ) x and setting H 2 k = 0, one 
obtains the expressions for all second order partial derivatives of s in terms of h and 5, which 
are analogous to (I14p . The consistency conditions are satisfied identically modulo (|14p . (|15p . 
Similarly, the condition H 3 k = specifies the right hand side of the fifth equation of the chain, 
etc. Although we know no direct way to demonstrate the non-obstructedness of this recursive 
procedure in general, there exists an alternative direct approach to the reconstruction of a 
chain from the function g(u l ,u 2 , u 3 ). To illustrate this procedure we consider a simple example 
g = u 3 — ^(u 1 ) 2 which automatically satisfies (|15p. The corresponding h, as specified by (|14p. 
is given by h = \i + au 1 + (3u 2 + 7U 3 + 5u 4 — ii 1 u 2 , which can be reduced to a canonical form 
h = u 4 — u l u 2 by redefining it 4 appropriately (this transformation freedom allows one to absorb 
arbitrary integration constants arising at each step of the construction). Thus, the first three 
equations of the chain are 

u\ = u 2 x , u\ = (u 3 - \{u l f) x , n 3 = (n 4 -nV) a .,.., (17) 

etc. Equations (|17p are nothing but the first three equations in the conservative representation 
of the Benney chain ([2]); notice that the variables u 1 in ([2]) and (|17p are not the same: u % in (fTj 
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are conserved quantities of the chain The recostruction of the remaining equations of the 
chain consists of three steps. 

(i) One introduces the corresponding PDE f)16[) : in our case this will be a potential form of 
the dispersionless KP equation, uu = u xy — ^u xx . It was demonstrated in [H] that the general 
PDE (|16p is integrable by the method of hydrodynamic reductions if and only if the function g 
satisfies the relations (I15p , 

(ii) One constructs a dispersionless Lax pair for the PDE (|16p . in our example it takes the form 

Pt = CjP 2 + Uxx)x, Py = {^P 3 + UxxP + Uxt)x\ 

the consistency conditions pt y = p y t are satisfied identically modulo the dispersionless KP equa- 
tion. The existence of such Lax pairs was established in [4U\ [T3] for any equation (|16p provided 
g satisfies the compatibility conditions (|15p . 

(iii) One looks for p as an expansion in the auxiliary parameter A, 

12 3 
U U U 



P ~ X A A 2 A3 •-' 
the substitution of this ansatz into the first equation of the Lax pair implies an infinite hydro- 
dynamic chain for the variables it*. The first three equations of this chain identically coincide 
with (|17p . The substitution into the second equation of the Lax pair produces a commuting 
chain (one has to set u 1 = u xx , u 2 = u x t). Both chains possess infinitely many hydrodynamic 
reductions since this is the case for the generating equation (|16p . Thus, the Haantjes tensor will 
automatically vanish (as demonstrated in Sect. 5). 

This procedure has been successfully implemented in [44J for various classes of solutions of 
the system ([15]) . namely, cases (fl~B"P and (fl~9j) of the classification presented below (according 
to [10], the case (fT9j) is reciprocally related to the case ([TBI and, therefore, does not require a 
special treatment). The work on (|20p and (|2ip is currently in progress 



2.1 Integration of the system ( 1151) 

To explicitly calculate g(u l ,u 2 ,u 3 ) we will follow [ID]. The main observation is that the first 
six equations in (fT5|) imply that the function I/53 is linear, I/53 = a + flu 1 + 7U 2 + 5u 3 . If 
5^0 then, up to a linear change of variables, one can assume that l/<?3 = u 3 . Similarly, if 
5 = 0, 7^0, one can set I/53 = u 2 . If (5 = 7 = 0, j3 7^ one has I/53 = vr. The last possibility 
is 1/(73 = 1- Thus, we have four cases to consider: 

3 / 1 0\ /I 2\ ^ 3 /12\ 1 3 /12\ 

g = u +p{u ,u ), g = — +p(u ,u ), g = + p{u ,u ), g = Inn + p(u , u ); 

here the function p{u l , u 2 ) can be recovered after the substitution into the remaining four equa- 
tions (fT5j) . In each of these cases the resulting equations for p{v},u 2 ) integrate explicitly, see 
[40J, leading to the four essentially different canonical forms: 

g = + -^{Au 2 + 2Bu 1 ) 2 + Ce~ Au \ (18) 

u 3 f I A \ . 2, 2 B 2 B 2 „ -A/u 1 

9 = -J + \-r ~ TTTv) ( u ) + T^ u ~ 77-^3 ~ Ce > ( 19 ) 
u 1 \ n 4(u 1 ) 2 / (u 1 ) 2 A(v}) 2 

g = t i + \i 1 {u 1 )(u 2 ) 2 , (20) 
g = lnu 3 -lno-(/u\u 2 ) - - frjiu^du 1 . (21) 
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Here 77 (it 1 ) is a solution to the Chazy equation [3], 

rj" + 27/77" - 3/ = 0, (22) 

and ^(u 1 , u 2 ) is an elliptic sigma function in the variable u 2 whose dependence on u 1 is governed 
by the Chazy equation (see Case IV below). Details of the derivation of canonical forms (|18p - 
(|2ip can be summarized as follows. 

Case I. Substituting the ansatz g = v? + p(u 1 ,u 2 ) into (fT5j) one arrives at the equations 

pm = 2(p\ 2 -PX1P22), 

P\\2 = P\22 = V222 = 0. 



The last three equations imply 



p = ^A(u 2 ) 2 + (Bu 1 + D)u 2 + q{u l ), 

and the substitution into the first equation results in the linear ODE q'" + Aq" = 2B 2 . Up to a 
transformation of the form it 3 — > u 3 + au 2 + [3u l + 7 this leads to (| 18|) . 

Case II. Substituting the ansatz g = u^/u 1 +p(u 1 ,u 2 ) into (fTBT) one arrives at the equations 

2 2 

p ni = 2{p\ 2 -PUP22) + T^rp ( pi +p 2) ~ ^T^ 11 + 2 V\P22 - 2p 2 pn), 

2 4 
P112 = - — P12, 



2 2 

P122 = 7— fTn fP22, P222 = 0. 

The last three equations imply 

( 1 A A / 2x2 /'■D B \ 2 B 2 , U 

and the substitution into the first equation results in the linear ODE (u 1 ) 3 g /// + n 1 (6 , u 1 — A)q" + 
(6u 1 -2A)q' = whose basis of solutions consists of 1 , 1 /u 1 and e A ^ ul . Up to a transformation 
of the form u 3 — > -u 3 + au 2 + /3m 1 + 7 this implies (|19|) . It was observed in [3D] that the cases 
I and II are reciprocally related: under the change from x,t to the new independent variables 
X, T defined as dX = u dx + u 2 dt, T = t, and the introduction of the new dependent variables 
U 1 = U 2 = — \, £7 3 = — ^t, etc, the chains from the case I transform to the chains from 

u u u 

the case II, and vice versa. On the level of the corresponding equations this means that the 

change of variables L rl = ^ r , U 2 = — P = ^-f p, G = &l g transforms the equations 

for p from the Case I to the equations for p from the Case II. Equivalently, (fTHj) goes to (fl~9|) . 
Case III. Substituting the ansatz g = u 3 /u 2 + p(u 1 ,u 2 ) into (|15p one arrives at the equations 

2 2 . 2 , 4 . . 

^ m = -mvlPi + 2 ^i2 -PHP22) + -n{PlP\2 -P2P11), 

2 2 
2 2 

P222 — / 9N9 P2 9^22- 
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The last three equations imply 

, B + Cu 1 1 . i,, 2 , 2 

p = a + — - — + -r]{u l ){u 2 y, 

u z b 

and the substitution into the first relation results in the Chazy equation (|22h for r\. Elimination 
of the constants A, B, C leads to (f20|) . 

Case IV. Substituting the ansatz g = lnu 3 + p(u 1 ,u 2 ) into (I15p one arrives at the equations 

pin = 2(pi2 -P11P22) + 2(pi + p|)pn - 4pip 2 pi2 + 2pip 22 , 
P112 = 4pipi2 - 2p 2 pn, 
P122 = 2pn + 2pip 22 , 
P222 = 4pi2 + 2p 2 p 2 2- 



The general solution of the fourth equation can be represented in the form 

p = — In a (u 1 , u 2 ) — — J rj{u l )du l ; 



here a solves the heat equation 4<ti = 022 and rj is a function of u . It is convenient to introduce 
the new variable v{v},u 2 ) by the formula v = — (In (7)22- Taking into account the heat equation 
for a one has 

1 1 2 in X 

f 1 = ^22 - 7^ + 2( m<T )2 t; 2- 

Rewritten in terms of v, the third equation for p implies 

U22 = 6v 2 - 4vr] - 4r/, (23) 
' = d/du 1 , the second equation is satisfied identically and the first takes the form 

Q 

v \ = 4u 3 - Av 2 r] - 8vr)' - -77". (24) 

3 

This shows that v is a shift of the Weierstrass elliptic function in the variable u 2 . Since v = 
— (In (7)22 ) the function a is the corresponding theta function. Notice that ([23]) can be obtained 
as a result of differentiation of (|24p by u 2 . Thus, we have two equations for v: 

v\ = v 2 - vt] - 7]' + -(lncr) 2 U2, 

i>2 = 4u 3 — 4t> 2 77 — 8^77' — -j/'. 

3 

The condition of their consistency leads to the Chazy equation for rj. 
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3 Classification of diagonalizable chains of the type ( 1121) 



In this section we sketch the classification of diagonalizable chains of the type (|12p . The condition 
Hj k = implies the following formulae for second order partial derivatives of h: 

, /&4513 , h 4 g 2 3 h 4 g 33 

tlU — , "24 — , ^34 — ; "44 — U, 

53 53 53 

, /2251 , 513 (h 3 - 52) + 52351 + 51253 + 533 /*i 

ni3 = 7 1 ! 

72 53 

, , , {h~92)hl , (^3 - /l)523 + /2513 + ^2533 

"23 — 522 — 7l2 H 7 1 ) 

72 53 

, /2253 533 (fl +52 - 2/l3) 

«33 = 2^23 7 1 (25) 

72 53 

, 722 (7i 2 ~ 2f l92 + g\ + jjgQ (fljj - /!)( 522 - 2/i 2 ) + / 2 (2 5l2 - /n) + 2 52 3^2 
^22 = -„ h 



/?. 
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53/2 53 

/225l(/l - 52) 52251 + /2511 + 513^2 + 523^1 ~ 2ffl/l2 
53/2 53 



, 5i/22 , 251512 + (fi - 52)511 + 2/11513 - 5i/n 

"11 = 7 — I • 

5372 53 

By calculating the consistency conditions for the above equations we obtain the expressions for 
all third order partial derivatives of 5 and /. 
Equations for 5: 



5333 — 


25 3 2 3 

53 


5133 — 


2513533 
53 


5233 — 


2523533 
53 


5113 = 


25i 2 3 
53 


5123 = 


2513523 
53 


5223 = 


25 2 2 3 

53 



5112 

-2 



/ll (-5335251 + 53(52351 + 52513 - 53512)) + 2/f(5i 3 ~ 5335ll) + /l/ll(5335l - 53513 

/25§ 

2/i ((-53351 + 53513)512 + 523(51513 - 535n) + 52(~5i3 + 5335n)) + 5 2 35i 

hoi 

522(~5i533 + 25351513 - 5 2 5n) + 5i3(52 + /25i) + 512(25335251 - 25253513 + 5 2 5i2) 

/25l 

gn (-53352 ~ /253351 + /253513) - 2523(5153512 + 52(51513 - 535ll)) 

/25| 

/12 (53351 + 53 (-251513 + 535ll)) 
/25s 

(5335i + 53(535n - 25i5i 3 ))/22 + /2(/n (53351 - 53513) - 2513(52513 - 253512)) 

/25l 

533(251512 - 52511) + 523(53511 - 251513) + /l (533511 - 513) 



5l 
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, (9192933 + (53513 - 51533) 7l + (53512 - 52351 ~ 52513)53) 
5122 — J22 



+ 



/25 3 2 

2 (5235l - 53352251 + 53522513 + /25l3 + (53351 ~ 53513 )/l2 ~ 725335ll) 

5| 



(((/l - 2g 2 )533 + 293523)/l + (52 + /25l)533 + (53522 - 2^2523 ~ 72913)93)722 

5222 — 7 — 2 

/253 

| 2 ~/l(/l2533 + 523 ~ 533522) + /l2(5 3352 - 53523) + 53523522 + 52(523 ~ 522533) 



5l 



| /2(533(/n - 4912) + 4923513) 



5§ 

Equations for f: 

fll (533(7l + 52 + 725l) - 2fif 2 53523 + 93 522 - 7253513) - 27l7ll(53352 - 53523) 



7m : 

+7l2 



725 3 2 

53/11 + 2 (5l92933 + 7l (-93391 + 93913) + 93(~9239l - 92913 + 93912)) 



729 3 2 

722(93391 - 29193913 + 93511) 

725| 

fnfwgi - 72(7i7n933 + 711(93923 - 93392) + 2712(93391 - 93913)) 



(27) 



'112 



'122 



7: 



222 



729! 

-7 2 2 7ll933 + 722(7l293 + 72(~9339l + 93913)) 
7291 

7 2 2 25 2 - 2 7 2 2 7i2533 + 72722((7i - 52)533 + 53523) 



/25§ 

We have verified that the system (j26H . (I27p is in involution. We claim that the functions 7 and 
9 contain all the necessary information about conservative chains (|12|) . In particular, equations 
(I25D allow one to reconstruct the function h and, hence, the right hand side of the third equation 
of the chain. Similarly, the requirement H? k = reconstructs the fourth equation (i.e., provides 
an involutive second order system for the right hand side of the fourth equation), etc. Although 
we have verified directly that our procedure is non-obstructed up to order 5, we could not 
establish this property in general. Possible approaches to this problem could be 

- developing a cohomological approach to the reconstruction procedure in the spirit of |8j ; 

- establishing a Hamiltonian formulation for general integrable chains of the type (|12|); 

- establishing a link of integrable chains of the type (|12f) to (2 + l)-dimensional integrable 
PDEs, as explained in Sect. 2. Although there is a general belief that one can construct 
a one-to-one correspondence between integrable hydrodynamic chains and integrable (2 + 1)- 
dimensional quasi-linear PDEs (that is, any integrable PDE can be 'decoupled' into a pair of 
commuting chains), there exists yet no rigorous proof of this statement. 

3.1 Integration of the equations (126ft . (E 



Notice that the first six equations in (J26J) are exactly the same as in Sect. 2. Thus, there are 
four essentially different cases to consider. 
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Case I: g = u 3 + p(u 1 ,u 2 ). Substituting this ansatz into ([2"o]) . (|2~7|) one arrives at the following 

relations. 

equations for p: 

Pl2(2pi2 - fll) + 2pn(/i2 -P22) 

Pill = 7 > 

J2 

/22PII /22P12 /22P22 

P112 — 7 1 P122 — 7 ; P222 — 7 > 

J2 J2 J2 



equations for /: 



, /ll(/l2 - P22) + 2/12P12 - /22P11 
'111 = 



2 



, /22/ll , /22/l2 , /22 

J 112 — 7 ) 7122 — 7 j J222 — ~~ 7 • 

J2 J2 J2 

The last three equations for / and the last three equations for p lead, up to elementary changes 
of variables, to the two possibilities. 
Subcase Ii: 

f = s{u 1 )e u \ p = q(u 1 )e u2 . 

The substitution of these expressions into the remaining equations for pm and fm leads to a 
system of coupled ODEs for s(n 1 ) and (/(u 1 ): 



q 



,„ 2{{q'f- qq ")- q 's" + 2 q "s' 



... s"s' + 2s' q' - qs" - q "s 

s = . 

s 

Setting q = —s', the second equation will be satisfied identically while the first implies a fourth 
order ODE s""s + 3(s") 2 — 4s' s'" = whose general solution is an elliptic sigma-function: 
s = cr(u 1 ), here (In a)" = — p, (p') 2 = 4p 3 — c (notice that 02 = 0, 93 = c). Thus, as a particular 
case we have 

f = a(u 1 )e u \ p = -a'(u 1 )e u \ 

Subcase I2: 

/ = (an 1 + b)u 2 + siu 1 ), p = ^A(u 2 ) 2 + Bu l u 2 + g^ 1 ), 

a, 6, A, I?=const. The substitution of these expressions into the remaining equations for pm and 
/in leads to linear ODEs for s(n 1 ) and ^(n 1 ): 



2£ 2 - Bs" + 2(o - A)i 

aw 1 + 
(a - A)s" + 2a£ 



s'" 



au 1 + b 



These equations are straightforward to solve. One needs to consider two different cases: a = 
0, b = 1 and a = 1, 6 = 0. If a = 0, 6=1 then, up to unessential integration constants, we 
have 

,i 

2A V " ' "A' 



s = ae~ Aul , q = ^(u 1 ) 2 - a^e'^ 1 + Pe~ 2Aul 
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The case a = 1, 6 = leads to 



A-r ; ' v ; ' * 2(^i-i) 2 v y i 

Case II: g = u^/u 1 +p(u 1 ,u 2 ). Substituting this ansatz into (j26|) . ((27]) one arrives at the 
following set of relations, 
equations for p: 

r,P2 + fx + hV\ ~ 4p2/i . (/l - P2X/11 - 4pi 2 ) + 4pi (/12 - P22) . 
Pin — ^ 7 7 n 2 1 T7 r 

Pi2(2pi2 ~ + 2pn(/i2 - P22) _ 

/2 U 1 

2 \ . /ii-4pi 2 2pi/ 22 /22P11 
P112 = (Ji - P2) H j 1 f~7 H 7 — , 

r, h 2 /, \ 1 /22(P2-/l) . /22P12 

P122 = 2— -j— ^ + -r(/i2 - P22) H 77 1 7 — , 



/22P22 , /22 . 

equations for /: 



V222 — — 7 1 — r 

h u 



fw fl f\1 1 ^22 1 

/111 = -7-(/l2 -P22 r) + tt( 2 ?2 - 2/i +pi 2 n ) - ——(2pi +pxxu ), 

, /ll/22 /l2 , /l2/22 /22 , /I2 

7112 — 7 ^ — T) /122 — 7 r> J" 222 ~~ ~T- 

J2 U 1 f 2 U 1 j 2 

The cases I and II are reciprocally related (we thank Maxim Pavlov for pointing out this 
equivalence): under the change from x,t to the new independent variables X,T defined as 
dX = v}dx + fdt, T = t, and the introduction of the new dependent variables U 1 = U 2 = 

— "77 , U 3 = — nr. etc, the chains from the case I transform to the chains from the case II, and 
vice versa. On the level of the corresponding equations this means that the change of variables 
U 1 = \, U 2 = — \, F = — -4-, P = \f — p transforms the equations for p, f from the Case I 
to the equations for P, F from the Case II. Thus, we will not discuss this case in any more detail 
here. 

Case III: g = v?jv? + p(u l ,u 2 ). Substituting this ansatz into (|26p . (|27p one arrives at the 
following set of relations, 
equations for p: 



2p 2 Pi(4pi2 - /n) + 4pn(/i - y> 2 ) Pi2(2pi2 - /11) + 2pn(/i2 -P22) 
/ " U " / 2 (n 2 ) 2+ n2/ 2 



/22P11 , n Pn 2/2P1 + u 2 / 22 (pi + n 2 pi 2 ) 



P112 = — 7 1- 2—5") P122 



/ 2 ' ~« 2 ' / 2 (n 2 ) 2 

J22(P2 ~ fl) ~ /2P22 . P22/22 . /l2 2 

P222 = 2 J- 1 h 2— - t^toC/i - P2); 

u 2 / 2 / 2 u 2 (V) 2 
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equations for /: 



/m = 7^(2(71 " P2) " P2 2 n 2 ) + -^(u 2 (f n + 2p l2 ) + 2 Pl ) - ^11 



/22/11 , /11 , /22/12 , A2 /: 



j 22 

112 — 7 1 5") 7122 — 7 ) 7222 — ~7 J ■ 

h U z f 2 h u z 

The last three equations for / and the last three equations for p lead to the two essentially 
different possibilities: 
Subcase IIIi : 

f = s{u 1 )(u 2 )\ p = q{u l ){u 2 ) k +\ 

fc=const. The substitution of these expressions into the remaining equations for pm and hn 
leads to the coupled system of ODEs for s(u 1 ) and (/(n 1 ): 

h 4-9 

III ' frtfl 1 n\/ '\2 r) / 7 I 1 \ // I II , r, II l\ 



ks 



-{2{k + 2){q'Y - 2{k + l)qq" - q's" + 2q"s' 



s'" = ^(s'V - k(k + 1)A + 2ks'q') -(k- l)q". 
ks 

Notice that under the substitution s = k = 1 the equation for q reduces to the Chazy equation 
(E2D for qiu 1 ) = ^(n 1 ). 
Subcase III2: 

/ = (an 1 + b) Inn 2 + ^(n 1 ), p = (/(n^n 2 + a(u 2 Inn 2 — n 2 ), 

a, 6=const. The substitution of these expressions into the remaining equations for pm and /m 
leads to the coupled system of ODEs for s(u 1 ) and (/(n 1 ): 



„, _ 8(q') 2 - 4qq" - 2q's" + Aq"s' 
au 1 + b 



Q 



_ W -q) + 2aq' 

1 1 J. T? • 

an 1 + b 

Case IV: g = Inn 3 + p(n 1 ,n 2 ). Substituting this ansatz into (|26p . (|27|) one obtains 
equations for p: 

p 22 (2p 2 - 2p n ) +Pii(P2{2p 2 ~ 4/i) + 2(/ 2 + / 12 + f 2Pl )) +p l2 (2p 12 - hi + 4pi(/i -p 2 )) 
Pm = 7 

72 

Pl(2pi/l2 + (/l ~P2)/ll) 
/2 

/22(Pll — Pi) . /„ , n . / , s 

Pll2 = 7 hpi(4pi2 - /ll) + 2pn(/i -P2), 

J2 

Pl22 = -7-(pi(/l -P2) + P12) + 2pi(p 2 2 - /l2) + 2/2P11, 
J2 

«/*22 2 
P222 = -7-(P22 + /l(2p2 - /l) - p 2 ) + 2(p 22 - /i 2 )(P2 - /l) - /22P1 + /2(4pi2 - hi); 
J2 
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equations for /: 



/in = 4r(/i +P% + hPi -P22 - 2/1P2 + /12) + 2^(pi(/i -P2) +P12) + -t^(Pi -P11), 

J2 J2 J2 

«/"22 fx 1 y*22 [ /*12 

/ll2 = 7 1" - PZ) + 2/l2Pl, /l22 = 7 1" /2/ll + /22P1, 

J2 72 

/222 = f ~f + 2f 2 fl2 + / 22 (P2-/l) 
J2 



Although this system is in involution and reduces to the corresponding Case 4 in Sect. 2 under 
the substitution / = u 2 , we were not able to integrate it in general. Let us just mention that 
the last three relations for / imply the Monge-Ampere equation /n/ 22 — /12 = a i ul )f2- This 
suggests a separable ansatz / = s(u 1 )r(n 2 ). A simple analysis leads to the two possibilities: 

(i) / = s(u 1 )(u 2 ) k , p = I ^s'(n 1 )(u 2 ) fc+1 -Inn 2 , where s satisfies the equation ss" — ^rr(s') 2 = 
0. Without any loss of generality one can take s = (u 1 ) 1 ~ k . 

(ii) / = s{u l )e u , p = —s'{u l )e u , where s satisfies the equation ss" — (s') 2 = 0. Up to a linear 
transformation of u 1 one has s = e u . 



4 Classification of diagonalizable Hamiltonain chains 

It was observed in [27\ [28] that the Benney chain ([2]) can be represented in the Hamiltonian 
form PD , 

U * \ dx dx^ ) 3u' 

where i? 2 - 7 = (i— l)nj+j_2, and /i = (n 3 + (n 1 ) 2 )/2 is the Hamiltonian density. Further integrable 
examples can be constructed by looking at Hamiltonian densities in the form h = u 3 u 2 ), 
and imposing the constraint = 0. This implies the relations 

Pill (2 + U 1 P22) = P11P22 ~ Pl2, PU2 = P\22 = V222 = 0, 

which, up to a natural equivalence h — > a/i + au 1 + 6n 2 + c, lead to the Hamiltonian densities 

h = u 3 + a^ 1 ) 2 + /3uV + 7 (u 2 ) 2 + ^(n 1 ) 3 (28) 

where the constants a, /3, 7,(5 satisfy a single relation /? 2 — 4a7 + 125 = 0. These densities were 
first introduced in [29], where it was shown that the corresponding Hamiltonian chains possess 
an infinity of conservation laws which Poisson commute and form a complete system. One can 
prove that all other components Hj k , i > 2, of the Haantjes tensor are identically zero. 

The aim of this section is to characterize all densities of the form h{u l , u 2 , u 3 ) such that the 
Haantjes tensor of the Hamiltonian chain (I13p is zero. As in the previous Sections, the conditions 
H\ = provide expressions for all of the ten third order partial derivatives of h, the simplest 
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six of them being 



5/l§ 3 5/113/133 5h 23 h 33 

/i333 — "777 j "133 — — T77 > ft 233 — — T77 > 

2/l3 2/13 2ll 3 

kl13 = 2h 3 ' (29) 

3/li3/l23 + 2h 33 h 12 
^123 = 



2fe 

'23 



, 3/^3 + 2/133/122 

«-223 " 



2/l 3 

the remaining expressions for hm, /1112, /ii22> ^222 are n °t written out explicitly due to their 
complexity (we have used Mathematica to manipulate with these expressions). It was verified 
that the system for h is in involution, and its solution space is 10-dimensional. In what follows, 
(|29p refers to the full system of ten equations. It was demonstrated recently in [21] that the 
relations (|29p imply the vanishing of all remaining components of the Haantjes tensor, that is, 

Hj k = H) k = for any i > 2. 

Before proceeding with a detailed analysis of the equations (|29|). let us make a digression on 
conservation laws of the chains (|13f) . First of all, any Hamiltonian chain (|13l) possesses three 
conservation laws of the form 

u\ = (u 1 h 2 + 2u 2 h 3 ) x , 
uf = (u 1 ^ + 2u 2 h 2 + 3u 3 h 3 - h) x , 

and 

ht = (u 1 /!!^ + 2u 2 hih 3 + u 2 ^! + 3u 3 h 2 h 3 + 2u 4 /i|) x , 

which correspond to the conservation of the Casimir, momentum and the Hamiltonian, respec- 
tively. Let us require the existence of an extra 'higher' conservation law of the form 

P(u l , u 2 , u 3 , u A ) t = Q{v},u 2 ,u 3 , u 4 , u 5 ) x (30) 

whose density P depends on the first four coordinates u 1 (the structure of the flux Q follows 
from the equations of the chain and does not constitute an additional restriction). We have 
obtained the following 



Proposition. The relations (2S\) are necessary and sufficient for the existence of an extra 
conservation law of the form \30\) . 

The proof is computational, see Appendix 2 for details. In fact, we have the following 
stronger result 

Theorem 4 The relations \29\) imply the existence of a generating function which gives rise 
to an infinity of conservation laws. To be more precise, we claim that for any n there exist n 
linearly independent conserved densities which are functions of the first n coordinates u l . All 
higher conservation laws are polynomial in n 4 , u 5 , u 6 , .... 

The proof of this statement is given in Sect. 4.2 below. 
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4.1 Integration of the system ( 1291) 

The first three equations i 
consider, s = and s ^ 0. 



The first three equations in f|29|) imply that /133 = sh^ 2 , s = const. Thus, there are 2 cases to 



Case I: s = 

In this case h is linear in u 3 and the equations (I29D imply 

7 ^ /I 2\ 

This ansatz can be simplified by utilizing the canonical transformations [31 



u 1 = An 1 , -u 2 = -u 2 , u 3 = — u 3 , u 4 = — j« 4 , (31) 

A A^ 

etc, A=const, and 

u 1 =u 1 , u 2 =u 2 + su 1 , u 3 = u 3 + 2su 2 + sV, (32) 



etc, s=const. Both transformations preserve the Poisson bracket specified by (j!3f) . Hence, they 
can be used to simplify the Hamiltonian. Suppose, for instance, that 6^0. Then, up to the 
second canonical transformation, one can assume the ansatz h = u 3 /(c + u 2 ) 2 + p{v} ,u 2 ). If 
6 = then h = u 3 / (c + u 1 ) 2 + p(u , u 2 ). Thus, there are two subcases: 

Subcase Ii: h = u 3 / '(c+u 1 ) 2 +p{u 1 , u 2 ) . The substitution of this ansatz into the four remaining 
equations for h implies the following system for p: 

16(u 2 ) 2 + (c + u 1 ) 3 (8(u 2 ) 2 p 2 2 + 4(3c + u l )u l pi 2 ) + 2(c + u r ) 4 (5c + v})pu 

PlU ~ c ( c + n l)5( 2 + n l( c + u l)2 p22 ) + 



(C + U l ){c 2 p\ 2 - c(c - 3u 1 )p22Pl 



c(2 + u 1 (c + u 1 ) 2 p2 2 ) 
2u 2 - c(c + u l ) 3 pi2 2 + c(c + u 1 ) 2 p22 

Pll2 = 2 r— ^ , P122 = 1— jTo , P222 = 0. 

c^c + u 1 ) 4 c(c + u L ) A 

The last three equations lead to p(v}, u 2 ) = ~^^rp - (^ 2 ) 2 + Jp^r + liu 1 ) where a and f3 are 
arbitrary constants. The substitution into the first equation gives a linear ODE for q, 



(ac 2 - Zv^ac - 4)g" - (c + u x )(l + cu l a)q"' 



c 2 (3 2 



2(c + u 1 ) 3 ' 

Without any loss of generality one has 

, t . m n 

) = 7—, — Tvl + 



(c + u 1 ) 2 c + u 1 

where the constants m, n satisfy a single relation 2(1 — c 2 a)n — Qcam + \c 2 (3 2 = 0. Ultimately, 
we have Hamiltonian densities 

o . . i, l + acic + u 1 ), n, 9 du 2 m n 

h = u 3 /(c + u 1 ) 2 + j 1 2 ' (u 2 ) 2 + + + — - -. 

c(c + u L ) A c + u L (c + u 1 )^ c + u L 
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Notice that the constant (5 can be eliminated by the second canonical transformation. 

Subcase l 2 : h = u 3 / (c+u 2 ) 2 +p{u 1 ,u 2 ). The substitution of this ansatz into the four remaining 
equations for h implies the following system for p, 

{u x ) 3 {c + u 2 )p ul - 2c{u 1 ) 2 p 11 - 2cu 1 (c - u 2 )p 12 - (5c - 3u 2 )(c - u 2 ) 2 p 22 



V222 

V\22 
V\\2 



(c — u 2 ) 3 {c + V?) 
(V ) 2 (c + u 2 )p lu - 2cu 1 p n - 2(2c 2 - 3cu 2 + (u 2 ) 2 )pi 2 

(c — u 2 ) 2 {c + u 2 ) 
u l {c + u 2 )pm + (u 2 - 3c)pn 
(c — u 2 ){c + u 2 ) 



along with one more relation for pm which we do not write out due to its complexity. The first 

three relations imply p(u l ,u 2 ) = ^—^ — ( c +jffi ~ " where g{f]) is an arbitrary function 

of a single variable 77 = . The result of the substitution of this ansatz into the remaining 
equation for pm factors into a product of two terms, leading to the cases (i) and (ii) below: 
(i) The function g satisfies a first order ODE 

4c 2 r] 2 g' - \2<?r)g - 1 - ar] + 2cr] 2 = 0, 

the solution of this equation is 51(77) = /i.77 3 + — — 16 * 2 . This results in the Hamiltonian 
densities of the form 

u 3 v}u 2 u 2 + c/3 (it 1 ) 3 1 (c — u 2 ) 3 ( u l a c - u 2 



(c + « 2 ) 2 (c + 7J 2 ) 2 (c + U 2 ) 2 (c + « 2 ) 2 4c(c + 7i 2 ) 2 \C — Ii 2 3c 4C7J 1 

(ii) The function 5 solves a third order ODE 

(4c 2 7 ? 2 y-12c 2 7 ?5 -l-a7 ? +2c7 ? 2 )/ / + (a-4c7 ? +12c 2 5 +4c 2 7 ? y-2c 2 r ? 2 /)/ + (4c-8c 2 c/ / )5'-^ = 0. 

Remarkably, this complicated equation trivializes after being differentiated by rj once, taking 
the form 

{Ac 2 r, 2 g' - 12c 2 rjg - 1 - ar) + 2cr] 2 )g"" = 0. 

Since the possibility when the coefficient at g"" equals zero was considered above in the case (i), 
we conclude that g must be a cubic polynomial, 

g = u + ur] + 7?/ 2 + 5r] 3 

where the constants satisfy a single relation 125 — 8cu + 16c 2 [y 2 — 37//) — 47a +1 = which can 
be obtained by substituting back into the original ODE. This leads to Hamiltonian densities of 
the form 



u 



3 



h = 7 ^7 + 



uV - era 2 - <f + n(c - u 2 ) 3 + vv}(c - u 2 ) 2 + 7(n x ) 2 (c - u 2 ) + 5{u v f 



{c + u 2 ) 2 {c + u 2 ) 2 

Notice that the particular case c = results (up to obvious equivalence transformations and 
relabeling of constants) in Hamiltonian densities of the form 



u 3 {u 1 ) 2 V 1 Ju 1 ) 3 
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where a, (3, 7, 5 satisfy a single relation j3 2 — 4q7 + 125 = 0. 

Remark. The apparent similarity of the cases (j28[) and (I33p is not accidental and manifests 
an important reciprocal invariance of the class of Hamiltonian chains (|13p . Recall that the 
conservation of momentum reads 

uf = {u l hi + 2u 2 h 2 + 3u 3 h 3 - h) x . 

Let us change from x, t to the new independent variables X, T where 

dX = u 2 dx + (u x hi + 2u 2 h 2 + 3u 3 h 3 - h)dt, T = t. 

It is known that reciprocal transformations of this type preserve Poisson brackets of the form 
(|13p . see e.g. [9] (it is crucial that u 2 is the momentum of the corresponding Poisson bracket). 
One can verify directly that performing the above change of independent variables and intro- 
ducing 

H =^ U ' = ^ U2 = ^ U3 = ^ t/B= (^' (34) 

one arrives at the system which takes the original form (|13l) in the variables X, T, U n , H. Thus, 
the above reciprocal transformation is canonical. One can verify directly that the change of 
variables indeed identifies ([28]) and (|33j) . 

Case II: s ^ 

In this case the elementary integration gives 

h = ^u 3 +p(u\ U 2 )) 1 / 3 + q(u\u 2 ), 



and the substitution into the last three equations (I29h implies that q is linear. Up to the 
equivalence h — > ah + au 1 + bu 2 + c we thus have 

h= {u 3 +p(u\u 2 )) 1 / 3 . 

The substitution of this ansatz into the remaining equations for h implies a complicated system 
of third order PDEs for piu 1 ^ 2 ). A useful observation is that this system is invariant under 
a 3-parameter group of point symmetries which is generated by the two canonical transforma- 
tions (|3ip . (|32p and the reciprocal transformation (I34p . The infinitesimal generators of these 
symmetries are 

Xi = u l d u i - pd p , X 2 = u x d u 2 - 2u 2 d p , X 3 = (u 1 ) 2 ^ + u l u 2 d u 2 + {2,pu l + 2(u 2 ) 2 )d p ; 

they satisfy the commutator relations 

[X U X 2 ] = X 2 , [X x ,Xz] = X 3 , [X 2 , X 3 ] = 0. 

These symmetries suggest a change of variables which considerably simplifies the equations for 
p. The idea is to choose new coordinates such that the symmetry generators assume the simplest 
possible form. Introducing £ = — rj = s = ^ , we have 

Xi = - rjd v - 4sd s , X 2 = d v , X 3 = d^. 
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Setting s = s(£, 77) we have p = s(£, ^(it 1 ) 3 — t . In terms of s(^, rj) the equations for p assume 
a remarkable symmetric form 

- s m s vv = 4s ? , (36) 

Srir/riS^ S^^S-q — Sq^S^, (37) 



1 



•^rpyS^S^ + S^-qS^S^ 4s^ — 5.^(125 + Sg^ + S^Sjjjj), (39) 

s «£ s w s ? + s &v s tt s v + 4s | = ^ s «( 12s + s |r? + s « s w); (40) 



in the process of derivation of these equations we have assumed that s^, and s^s 2 — s vv s 2 are 
nonzero: these expressions appear as denominators in the intermediate calculations. Particular 
cases when some of these expressions vanish will be discussed below. 

To solve the equations (f35|) - (jlOj) we proceed as follows. Differentiating (|3T|) . (j38|) and using 
([35]) . ([36]) we obtain four relations among the fourth order partial derivatives of s, 

^Sijtjijt) + SfjSj:^^ = 0, S^S^rjrir) ^17(^^7777 4) = 0, 

s ?( s ^w - 4) + Sr?s^ = 0, S£S^ V + s,,s^ = 0, 

which can be parametrized as 

s V vvv = <li s tvrm = <l r i s ttrm - 4 = qr 2 , s^ v = qr 3 , = qr 4 , (41) 

r = —s^/s v . The further analysis leads to the two possibilities. 
Subcase Hi : q = 0. In this case s(£, 77) is a polynomial of the form 

s = £V + < 3 + bfv + c£f? 2 + drf + af + + jr] 2 + /if + wq + e. 

The substitution into the remaining equations for s implies the following relations among the 
coefficients: 

(3 = be — 9ad, fi = ac — 3a7, v = 67 — 3ad, 12e + f3 2 + 4aj = 4(bu + c/j). 

Notice that the corresponding p(v}, u 2 ) = s{—\, ^t)(m 1 ) 3 — ^1 will be a cubic polynomial in 
u 1 ,^. A particular example from this class with the Hamiltonian density h = (u 3 + r) 1 / 3 was 
discussed in [29J. It corresponds to the case where a = — r and all other coefficients of s{^,rj) 
are zero. 

Subcase II2: q ^ 0. Then the consistency conditions of (]4ip imply the relations q^ = (qr) v , r^ = 
rr v . Taking into account that r = —s^/s v we have s^s 2 — = 0. This case is discussed 

below. 

There are two more possibilities one needs to consider to complete the classification (notice 
that the equations (f35|) - (|4"0l) can no longer be used since they were derived under the assumption 
that certain expressions do not vanish). 

Subcase II3 : = or s v = 0. A simple analysis leads to Hamiltonian densities of the form 



h = [u 3 + a(u) 



us I" - 



2X2 \ 1/3 



It 1 
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a =const. 

Subcase II4: s^s 2 — s w s| = 0. Setting r = — one obtains 

r% = rr v , s$ = -rs v . 

Calculating partial derivatives of s with the help of the above relations and substituting them 
into the conditions Hj k = we obtain r w = 0; without any loss of generality one can set r = — |. 

This implies — r/s^ = 0, therefore, s = s(z), z = £77 = — u 2 /{u 1 ) 2 . For s(z) we obtain the 
ODE 

8zVV + 8zs" S ' - 4z 2 (s") 2 - (s') 2 - 12s = 
which linearizes after being differentiated by z once. Its general solution is given by the formula 

s{z) = z 2 + a + (3(-z) + l(-z) 1 ' 2 + 5(-zf/ 2 

where the constants a, (3, 7, 6 satisfy a single quadratic relation 12a + (3 2 — 3^/5 = 0. This results 
in Hamiltonian densities of the form 



h = (u 3 + a(n 1 ) 3 + (5u l u 2 + 7 (n 1 ) 2 (n 2 ) 1 / 2 + 5(u 2 f' 2 ) 



1/3 



4.2 Generating functions for conservation laws 

In this section we demonstrate that any Hamiltonian chain (|13p whose density h satisfies the 
system (I29p . possesses a generating function which provides an infinity of conservation laws. To 
illustrate the method of generating functions let us consider the Benney chain (|2|) and introduce 
the generating function 

12s 

u u u , , . 

\ = p + — + ^ + + ... (42) 
p p z p^ 

which was shown in [17] to satisfy, by virtue of the fundamental relation 

X t - p\ x = X p [pt - {p 2 /2 + u l ) x ] . (43) 
This relation provides an infinity of conserved densities H n (u) defined by the equation 

p t = (p 2 /2 + u 1 ) x 

where one has to substitute the expression for p(X) obtained from pSjk p = X — — "fir - 
Explicitly, one gets Hi = u 1 , H 2 = u 2 , H 3 = u 3 + (n 1 ) 2 , etc. 

Given an arbitrary Hamiltonian density h(u l , u 2 , n 3 ), a generating function is sought in the 
form 

00 f, 

X = ^{p,u\u 2 ) + Y^; (44) 
k=i P 

notice that the dependence on higher moments n 3 , u 4 , u 5 , etc, is exactly the same as in (|42p . As 
demonstrated in [44J, the generating function (|44jl has to satisfy the Gibbons- type relation 

A t - (2ph 3 + h 2 )X x = X p [pt - (p 2 h 3 + ph 2 + h x ) x ] , (45) 



22 



which reduces to ([32]) for the Benney Hamiltonian h = (n 3 + (n 1 ) 2 )/2. Substituting ([33]) into (|45|) 
one obtains, by virtue of (|13H . the following relations among the first order partial derivatives 

1pl,1p2,1pp- 

3u 3 

(V>j> r)(p 2 hi3 +phu + /in) + ipi(2u 2 hi 3 + u l h l2 - 2ph 3 ) + V 2 (3u 3 /i 13 + 

1 3u 3 
2u 2 hu + u x h n ) + ^{3u 3 h X2 + {2u 2 + )h n ) = 0, 

3u 3 

(.4>p ~ —r)(p 2 h 2 3 +ph 22 + h 12 ) + M2u 2 h 2 3 + u l h 22 + 2h 3 ) + 4> 2 {3u 3 h 23 + (46) 
pi 

1 3u 3 

2u 2 h 22 + u x h X2 - 2ph 3 ) + -^{3u 3 h 22 + {2u 2 + )h l2 ) = 0, 

pi p 

3v? 

rXAss + ph 23 + /lis) + M2u 2 h 33 + u l h 23 ) + ip 2 (3u 3 h 33 + 2u 2 h 23 + 

pi 

1 3u 3 

u l h 13 + 2h 3 ) + -^(3u 3 h 23 + (2u 2 + )h l3 - 2ph 3 ) = 0. 

p^, p 

These relations uniquely specify the partial derivatives ipi,ip 2 ,i/j p : 

ipi = 

-— [4p 2 h 3 (ph 3 h 33 + h 23 (2h 3 + v}h 13 ) - /i 33 (V/ii 2 + u 2 h 22 )) 

p(2h 3 (3u 3 (h 33 h 22 - h 2 23 ) + 2h 3 {h 22 + h l3 ) + u\h 13 {2h 22 + h 13 ) - 2h 12 h 23 - h n h 33 )) + 
(n 1 ) 2 (/ii 2 (/ii2 + h 33 - h l3 h 23 ) + h 23 (h u h 23 - hi 2 h 13 ) + h 22 (h\ 3 - h\ih 33 ))) 
h 3 (6u 3 (h 13 h 23 - h 12 h 33 ) - 2h 3 h 12 + v}(hiih 23 - h 12 h 13 )) + 2u 2 u 1 (h 12 (h 12 h 33 - h 13 h 23 ) + 
h 23 (h u h 23 - h 12 h 13 ) + h 22 (h\ 2 - /in/133))] , 



1P2 = 

-— [iphlh 33 p 2 + 4h 3 {h 23 (h 3 + u 2 h 23 + v}h 13 ) - h 33 (u 2 h 22 + u 1 hi 2 ))p+ 

(n 1 ) 2 (/ii 2 (/ii 2 /j33 - h l3 h 23 ) + h 23 (huh 23 - h 12 hi 3 ) + h 22 (h\ 3 - h u h 33 )) + 
2h 3 (h 13 (2h 3 + 2u 2 h 23 + u\h 13 + h 22 )) - h 12 (2u 2 h 33 + u l h 23 ) - u^u/tsa)] , 



% = 

-— [4hlh 33 u x p 2 + 4h 3 p(2h 3 (u 2 h 33 + u% 3 ) + uV(/^3 - /i 22 /i 33 ) - (v, 1 ) 2 (hwhva - h 12 h 33 ))+ 

2h 3 (2h 3 (3u 3 h 33 + 4u 2 h 23 + u 1 {2h 13 + M + (u 1 ) 3 {hi 2 {h 12 h 33 - h l3 h 23 ) + h 23 (h 23 h n - h l2 h l3 ) + 
h 22 {h\ 3 - huh 33 ))) + 4(u 2 ) 2 (/i 2 . 3 - h 22 h 33 ) + 4uV(/n 3 /i 23 - h 12 h 33 ) - 3«V(/i| 3 - /i 22 /i 33 ) + 
{u l ) 2 {h l3 {2h 22 + h X3 ) - 2h 23 h 12 - h u h 33 ))] . 
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A = 

4p 4 hjh 33 + Ap 3 h 2 (2h 3 h 23 + u 2 (h 2 23 - h 22 h 33 ) + u l (h 13 h 23 - h 12 h 33 )) + {{^fp 2 + 2u 2 u 1 p) 
(2h 3 (h 13 (hy 3 + 2/t 2 2) - 2h 12 h 23 - h n h 33 ) + (h 12 (h 12 h 33 - h V3 h 23 ) + h 23 (hnh 23 - h 12 h 13 ) + 
M^ia - h n h 33 ))) + 2/i 3 p 2 (3u 3 (/i22/i33 - ^23) + 2/13(2/113 + M) + 4h 3P {3u 3 {h 13 h 23 - h 12 h 33 ) 
2h 3 h 12 ) + 2h n {h 3 (2h 3 + u 1 /^ + 4u 2 /i 23 ) + 2(u 2 ) 2 {h 2 23 - h 22 h 33 )) - 
3u 3 ((/i 2 3 - h u h 33 )(2h 3 + u l h 22 ) + u 1 h 23 (hi 2 h 33 - h X3 h 23 ) + u^h^ihuh^ - h 12 h 13 )). 

These equations are consistent provided the density h satisfies the system ()29p . This proves 
Theorem 4: indeed, an infinite sequence of conservation laws results from the equation 

p t = (p 2 h 3 +ph 2 + h\) x 

where one has to substitute the expansion for p in terms of A obtained from (I44p . 

The calculation of ipipjU 1 ,u 2 ) can be summarized as follows: 
- one first integrates the equation for ip p , which appears to be rational in p, with respect to p. 
This defines ip up to a function of u , u 2 . 

— one fixes this function by a substitution into the equations for ip\,ip 2 . 

This procedure leads to the following explicit formulae for generating functions for the Hamil- 
tonian densities obtained in Sect. 4.1. 

Linear case. This is the simplest generalization of the Benney Hamiltonian with the density 

h(u\ u 2 , u 3 ) = u 3 + aiu 1 ) 2 + Qu l u 2 + 7 (u 2 ) 2 + ^(u 1 ) 3 , 

where the constants satisfy the condition (3 2 — Aa^f + 125 = 0. Even in this case the function ip 
is quite nontrivial, 



tp(p, u x ,u 2 ) 



(4a 7 - B 2 



arctan 



P + 2jp 



(l + 7-u 1 )(4a7-/3 2 )2 



The generating function is obtained by substituting this expression into (I44h . 



Subcase I\. In this case the Hamiltonian density takes the form 

■ 3 l + acfc + u 1 ). 9 , 9 Bu 2 

(u z Y + 



h(v},u 2 ,u 3 ) 



u 



(c + u 1 ) 2 



+ 



c(c + u 1 ) 



1\2 



+ 



m 



+ 



n 



c + u 1 (c + u 1 ) 2 c + u 



where the constants satisfy the single relation, 2(1— c 2 a)n— 6cam+^c 2 B 2 = 0. The corresponding 
function ip is given by 



3 

C2 



(3m) 



tanh 



2(c 2 a - l)u 2 + c(2(c 2 « - l)p - pc) - cfiu 1 
2(3cm)3(l + cau 1 ) 



Subcase I2 (i) - For the Hamiltonian density 



u 3 + u 1 u 2 -a(u 2 + §) + n(u 1 ) 3 (c-u 2 ) 3 ( u 1 



h(u,U,U 3 ) = ^rf h 



c + u 2 ) 2 4c(c + u 2 ) 2 \c — u 2 3c 4cm 1 



a c — u 
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the function ip is given by 



(3/z)! 



arctan 



(u 2 ) 2 - 2cu 2 + Acpu 1 + c 2 



4c(3/i)5(u 



1^2 



Subcase II3: For the Hamiltonian density 

/ l (n 1 ,n 2 ,n 3 ) = (n 3 + a(n 1 ) 3 -^)t, 



the function ^ is given by 

ij>{j>,u x ,u 2 } 



(3a 



— arctan 



1 2 
pu — u 



(3a) 2 (u 



l^2 



Other examples obtained in Sect. 4.1 lead to more complicated expressions for ip. 

5 Hydrodynamic reductions and the diagonalizability 

To illustrate the method of hydrodynamic reductions we consider the Benney chain ([2]), 

U t = U x + ^ u xj 

u t = U t + 2« 2 ^, 
u 4 = itjj. + 3u 3 ul,, 

etc. Following the approach of [19j [20] let us seek solutions in the form u % = u l (R \ . . . ,R m ) 
where the Riemann invariants R 1 , . . . , R m solve a diagonal system 

R\ = \ % {R)R l x . 

Substituting this ansatz into the Benney equations and equating to zero coefficients at R % x we 
arrive at the following relations: 

(47) 
(48) 
(49) 
(50) 

etc. Here u = u , <9j = dj^i, i = 1, ...,m (no summation!) The consistency conditions of the first 
three relations (f4T|) - (f4*9"j) imply 

didju = — —diU + — —djU, 



V 2 


= \%u, 




u 3 


= ((A 4 ) 2 


- u)diu, 


u 4 


= ((A*) 3 


- uX 1 - 2u 2 )diU, 


u 5 


= ((A 4 ) 4 


- u(A*) 2 - 2u 2 y - 3u 3 )d iU 



\i -A* 



A 1 - \3 



dj\ l diU + di\ J djU = 0, 
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X t djX l d i u + \ j di\ j djU + diudju = 0, 
respectively. Solving these equations for dj\ l we arrive at the Gibbons-Tsarev system 

J - A* J (A* - \i) 2 v ; 



It is a truly remarkable fact that all other consistency conditions (e.g., of the relation (|50p . 
etc), are satisfied identically modulo (|5ip . Moreover, the semi-Hamiltonian property Q is also 
automatically satisfied. Thus, the system (|5ip governs m-component reductions of the Benney 
chain. Up to reparametrizations R l — > p(R l ) these reductions depend on m arbitrary functions 
of a single variable. Solutions arising within this approach are known as multiple waves, or 
nonlinear interactions of planar simple waves. 

The above approach clearly applies to any hydrodynamic chain from the class C. Let us 
restrict, for instance, to the chains of the type ([6]). Looking for solutions in the form u l = 
u % (R , . . . , R m ) and substituting this ansatz into ([6]) we arrive at an infinite system of relations 
similar to ()47p - (l50|) . The first three of them imply the 'generalized Gibbons-Tsarev system' of 
the form 

djX 1 = (...)djU, didju = (...)diudjU, 

u = u , where dots denote complicated expression which are rational in A* with the coefficients 
depending on the chain under study (that is, on g, h, etc). Requiring that all other consistency 
conditions, as well as the semi-Hamiltonian property, are satisfied identically we obtain con- 
straints for the matrix V(u), see [lOj EU E2 EH EH EE] where a similar approach was applied 
to the classification of integrable multi-dimensional quasilinear systems. The main result of this 
section is the proof of the following theorem formulated in the introduction: 

Theorem 2 The vanishing of the Haantjes tensor is a necessary condition for the existence 
of infinitely many hydrodynamic reductions and, thus, for the integrability of a hydrodynamic 
chain. 

We will give two different proofs of this statement. Based on essentially different ideas, they 
may be of interest in their own right. 

First Proof: 

This proof is computational. Writing down the equations of the chain in the form u™ = V™u™ 
and substituting the ansatz u l = u % (R l , R m ) we arrive at an infinite set of relations 



v™d iU n = A*a 



u 



we point out that all summations here and below involve finitely many nonzero terms. Applying 
the operator dj, j / i, we obtain 

V™ k diU n d jU k + V™didju n = djXdiu" 1 + X^idju" 1 . (52) 

Interchanging the indices i and j and subtracting the results we arrive at the expression for 
didjU m in the form 

o , j P) \j V m V m 



A; A' ' ' A' A' A' A> 

Substituting this back into (f52l) we arrive at a simple relation 

N"tdiU n d<jU k 



d^dm 171 + di\ j d iU m = nk x : J 
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where iV is the Nijenhuis tensor of V. This can be rewritten in the invariant form 

dj mu + dl ^u= N ^^ (53) 
which implies the following four relations: 

(X^d 3 X%u + (X^XWjn = V2N{ X ^ U \ 

x i xw j x%u+ xWdiX'dju = N <y^™ iU) ■ 

For instance, the first relation can be obtained by applying the operator V 2 to (|53p and using 
VdiU = A^u. Thus, 

V 2 N(diU, dju) - VN(VdiU, dju) - VN(diU, VdjU) + N(VdiU, Vdju) = 0. 

The last relation can be rewritten in the form H(diU, dju) = where H is the Haantjes tensor, 
indeed, a coordinate-free form of the relation (llli is 

H(X, Y) = V 2 N(X, Y) - VN(VX, Y) - VN(X, VY) + N(VX, VY) 

where X, Y are arbitrary vector fields. Keeping in mind that (a) <9jU and dju are eigenvectors of 
the matrix V corresponding to the eigenvalues X 1 and X J , and (b) A* and X 1 can take arbitrary 
values, we conclude that H(X, Y) = for any two formal eigenvectors of the matrix V. 

Assuming that formal eigenvectors of V span the space of dependent variables u (this is 
true for all examples discussed in this paper), we obtain H = 0. In more detail, let X(X) = 
£ 2 M> C 3 (^)i •••)* be a formal eigenvector of the matrix V corresponding to the eigenvalue 
A. Let us assume that these eigenvectors span the space of dependent variables, that is, that there 
exist no non-trivial relations of the form Cj£ l (A) = with finitely many nonzero A-independent 
coefficients Cj. In other words, £ J (A) are linearly independent as polynomials in A. 

The condition H(X(X), X(/j,)) = 0, written in components, takes the form Hj k ^ (X)^ k (fi) = 
0; recall that, since the matrix V belongs to the class C, these sums contain finitely many terms 
for any fixed value of the upper index i. Taking into account the linear independence of £ J (A) 
and £ fc (^)) one readily arrives at H l - k = 0. As an illustration, let V be the matrix corresponding 
to the Benney chain ©. Then X(X) = (1, A, X 2 -u 1 , A 3 - n 1 A - 2u 2 , ...)*, so that each f*(A) is 
a polynomial in A of the degree i — 1. These eigenvectors span the space of dependent variables 
since polynomials £ l (A) are manifestly linearly independent. 

Notice that we have proved a more general result, namely, that the existence of sufficiently 
many two-component reductions already implies the vanishing of the Haantjes tensor H. Indeed, 
nothing changes in the proof if we set i = 1, j = 2 in the formula (|53|) , As demonstrated in the 
Theorem 3 below, one can further strengthen the result by reversing the above proof under the 
additional assumption of the simplicity of the spectrum of the matrix V. This will require the 
relations 

d X i d u _ N ( Vd ^ d,") - VN(d t u, dju) j _ Njydpi, d lU ) - VNjdjU, d.u) 

1 3 {X-Xi) 2 ' 3 1 (A^-AJ) 2 
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which can be obtained by applying V to both sides of ([55]) and solving for diXWj\i and dj\ l diU. 

Second Proof: 

Our first remark is that for the chains from the class C one needs to know only finitely many 
rows of the matrix V(u) to calculate each particular component of the Haantjes tensor. Let us 
fix the values of indices i,j, k and denote by C(i,j, k) the maximal number of rows needed to 
calculate H l - k (counting from the first row). We need to show that H l - k = 0. Let us consider 
an m-component diagonal reduction u l (R}, . . . , R m ), i = 1,2, .... Choosing the first m variables 
u , u m as independent, we can represent the reduction explicitly as 

u m+l =u m+1 (u\...,u m ), u m+2 =u m+2 (u 1 ,...,u m ),..., 

etc. Substituting these expressions into the first m equations of the chain we obtain an m- 
component system S m for u , . . . , u m , while the remaining equations will be satisfied identically 
(by the definition of a reduction). Notice that the Haantjes tensor of the reduced system S m is 
identically zero since the reduction is diagonalizable. Let us now choose the number m sufficiently 
large so that the first C(i,j,k) equations of the chain do not contain variables u m+1 ,u m+2 , . . . 
(one can always do so since any equation of the chain depends on finitely many u's, and m 
can be arbitrarily large). Then the first C(i,j,k) equations of the reduced system S m will be 
identical to the first C(i,j,k) equations of the original infinite chain. Hence, the corresponding 
components H l - k for the reduced system and for the infinite chain will also coincide. This proves 
that all components of the Haantjes tensor of the chain must be zero. 

A straightforward modification of the second proof allows one to show that the existence of 
an infinity of semi-Hamiltonian reductions implies the vanishing of the tensor P. This establishes 
the necessity of the conjecture formulated in the Introduction. 

We emphasize that the condition of diagonalizability alone is not sufficient for the integra- 
bility in general. This can be seen as follows. 

Example. Let us consider the chain 

u] = u 2 x +p(u l )ul, 

u t = ul+p{u l )u 2 x + u l u l x , 

u 3 t =u x +p(u 1 )u x + 2u 2 u 1 x , 

etc, which is obtained from the Benney chain Uf = V(u)u x by the transformation V — > V + 
p(u 1 )E where E is an infinite identity matrix and p is a function of nr. One can verify that 
the corresponding Haantjes tensor is zero (which is not at all surprising since the addition of a 
multiple of the identity does not effect the diagonalizability). A simple calculation shows that 
hydrodynamic reductions of this chain are governed by exactly the same equations as in the 
Benney case, the only difference is that now the Riemann invariants K 1 solve the equations 

Ri = (X i (R)+p(u 1 ))Rl 

The semi-Hamiltonian property is satisfied if and only if p" = 0. Thus, we have constructed 
examples which possess infinitely many diagonal hydrodynamic reductions none of which are 
semi-Hamiltonian (if p" ^ 0) . 

Theorem 2 can be strengthened under the additional assumption of the simplicity of the 
spectrum of V. 
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Theorem 3 In the simple spectrum case the vanishing of the Haantjes tensor H is necessary and 
sufficient for the existence of two-component reductions parametrized by two arbitrary functions 
of a single argument. 

Proof: 

The necessity part is contained in the first proof of Theorem 2. To establish the sufficiency one 
has to show that the vanishing of the Haantjes tensor implies the solvability of the equations 

Fdiu = A^iu, Vd 2 u = X 2 d 2 u, (54) 
d!d 2 u m = ^^diu™ + JTT^ 9 ^ + n ^_^ n dm n d 2 u k (55) 



and 



o , 2 fl „ _ NjVdiu, d 2 u)-VN(d lU , d 2 u) 

o\A o 2 u — (x^-X 1 )' 2 > 

o ilQ „ _ NTVchu, diu)-VN(d 2 u, giu) 
2 A 0\U — 1 



(56) 



(A!-A 2 ) 2 

which govern two-component reductions. Our first observation is that the vanishing of the 
Haantjes tensor implies that the vectors N(Vdiu, d 2 u) — VN(c\u, d 2 u) and N(Vd 2 u, d±u) — 
VN(d 2 u, diu) are the eigenvectors of V with the eigenvalues A 2 and A 1 , respectively. By the 
assumption of the simplicity of the spectrum, they are proportional to c^u and <9iu. Thus, 
equations (|56|) reduce to a pair of first order PDEs for A 1 and A 2 , 

dlX2= (\i-\*)2> &2>} = (Ai- 2 A 2 ) 2 ' (57) 

here k\ and k 2 are the corresponding coefficients of proportionality. The relations (I54|) allow one 
to reconstruct all components ii 2 ,u 3 ,... of the infinite vector u in terms of its first component 
u 1 . Finally, the equations (|55|). which are nothing but the consistency conditions of (|54p . reduce 
to a single second order PDE for the first component u 1 of the infinite vector u. Thus, relations 
dMI) - dSHJ) reduce to a pair of first order equations ([57]) plus one second order PDE for u . 
Up to reparametrizations R 1 —> f 1 (R 1 ), R 2 — > f 2 (R 2 ), their general solution depends on two 
arbitrary functions of a single variable. 



6 Conclusion 

We have proposed a simple and easy-to- verify necessary condition for the integrability of hydro- 
dynamic chains based on the vanishing of the Haantjes tensor. Conservative and Hamiltonian 
examples are discussed, illustrating the general approach. We conjecture that all examples aris- 
ing in Sect. 3 and 4 from the requirement of the vanishing of the Haantjes tensor (and, thus, 
satisfying a necessary condition for the integrability) also satisfy the following properties: 

(i) they possess infinitely many m-component hydrodynamic reductions parametrized by m 
arbitrary functions of a single variable (thus, they are integrable in the sense of Definition 3); 

(ii) they belong to infinite hierarchies of commuting hydrodynamic chains. 

Our method leads to an abundance of new examples of hydrodynamic chains which require 
a further detailed investigation. We will address these issues in the future. 

We hope that the results of this paper will find important applications in the theory of 
infinite-dimensional Frobenius manifolds (yet to be constructed). 
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Appendix 1: the invariant formulation of the semi-Hamiltonian 
property 

For m-component systems (|7|) there exists a tensor object which is responsible for the semi- 
Hamiltonian property. First of all one computes (1, 3)-tensors M and K, 

M 8 kij = + N^T k N* - N; q N? k v] - N; q N p kj vf - N% p N? q v] - N s kp N^i 

and 

KUj = K d *K ~ %^ N 'j + N?d uP bi - N* p d u ^ + N* kp d uj q 
+Kd u ,Nj k - b^N^ + Nf^bf - N? p d u3 bl + Nf p d u ^ 

+K d ^ N ki ~ ^ N ki + K 9 ^ - N! p d u >% + N! p d tti %; 
here b = v 2 , that is, 6* = v p v v - . Using M and K one defines a (1, 3)-tensor Q as 

Qkj = < K UA + < K ^I - <vl KS pij - K i m <v) 

+4v p k M^ - 2M^ pj vf - 2Mi %p v). 
Ultimately, one introduces a tensor P, 

PL, = v pQ P kqj v i + KQlirf - v^Qlij - Qtp^v]. (58) 

Theorem 5 fJJjj A diagonalizable system (f7|) is semi-Hamiltonian if and only if the tensor P 
vanishes identically. 

These objects can be calculated using computer algebra. Notice that they are well-defined for 
hydrodynamic chains from the class C: all tensor operations will involve finite summations only. 
Invariant coordinate- free definitions of the above tensors can be found in [41] . 

It was pointed out in [46j (Proposition 5, p. 19) that any strictly hyperbolic conservative 
system 

ui = f(u) x , i = l,...,m (59) 

with the zero Haantjes tensor is necessarily semi-Hamiltonian. The shortest proof of this state- 
ment known to us can be summarized as follows. Let us first rewrite (|59|) in the Riemann 
invariants, R\ = X l (R)R x , where A* ^ X 3 due to the strict hyperbolicity. The conserved densi- 
ties u(R) satisfy an over-determined system of second order linear PDEs |51| . 

d A* d-X 3 
did 3 u = \A U + A i l _ A j dj M ' * + 3- (6°) 

The consistency conditions dk(didju) = dj(didku), i ^ j ^ k, lead to linear relations among 
the first order derivatives of u of the form ^ c n d n u = where c n are certain functions of A's. 
It is well-known that the vanishing of all coefficients c n , that is, the involutivity of the linear 
system (I60p . is equivalent to the semi-Hamiltonian property ([9]). In this case we have infinitely 
many conserved densities parametrized by m arbitrary functions of a single variable. It turns 
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out that the requirement of the existence of m functionally independent solutions of the linear 
system (|60l) is already sufficient to conclude that all coefficients c n must be zero. Indeed, if a 
relation of the form ^ c n d n u = is satisfied by m functionally independent it's, the vector c n 
will be a zero eigenvector of the corresponding Jacobian matrix which is non-degenerate. Thus, 
c n = and we have shown that 

(m conservation laws) + (Riemann invariants) =>■ (oo of conservation laws). 

It would be of interest to obtain a direct tensor proof of the above result by showing that 
the vanishing of the Haantjes tensor H for a conservative system (I59D implies the vanishing of 
the tensor P. Such a proof would then generalize to hydrodynamic chains, for which we have 
yet no rigorous definition of Riemann invariants. 

Appendix 2: higher conservation laws for Hamiltonian chains 

Let us assume that there exists a conservation law of the form 

P(u 1 ,u 2 ,u 3 ,u A ) t = Q(u 1 ,u 2 ,u 3 ,u 4 ,u 5 ) x . 
Using (Tl3|) and collecting coefficients at u x we obtain the expressions for Qj = dQ / d u i , 

2 3 

Qi = 5u 5 P 4 h 13 + J2( 4uip S-ih,(i+i) + ^Ps-ihij) + J2( 3u3p 5-3h ld + 2u 2 P 4 - j h hj ) + 
i=i j=i 

Pih 2 , 

2 3 

Q 2 = 5u 3 P 4 h 23 + ^(4u 4 P 5 -^ 2 ,(i+i) + u 1 P 3 ^h 2 , i ) + J2( 3u3p 5-j h 2,j + 2u 2 P^ j h 2J ) + 
i=i j=i 
P 2 h 2 + 2P 1 h 3 , 

2 3 

Q 3 = 5u 5 P 4 h 33 + ^(4u 4 P 5 -i/i 3 ,(i+i) + n 1 P 3 - l /i3, 4 ) + ^(3u 3 P 5 -i^3,j + 2u 2 P^ j h 3j ) + 
i=i j=i 
P 3 h 2 + 2P 2 h 3 , 

Qi = h 2 P 4 + 2h 3 P 3 , 

Q 5 = 2/13P4. 

The consistency conditions Qij = Qji imply the following expressions for second order partial 
derivatives of P, 
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6u 4 P 4 /i 2 3 + 9u 3 P 4 /ii 2 /ii3 + 3u 2 P 4 (/i 2 2 + 2/113/111) + 3^4/111/112 + 2P 3 h 3 h u 

12n 4 P 4 /i2 3 /ii3 + 9-u 3 P 4 (/i 22 /ii3 + h 23 h 12 ) + 6u 2 P 4 (h 12 (h 22 + h 13 ) + h 23 h u ) 

Pl2 ~ 4^1 + 

3tt 1 P 4 (/tf 2 + h 22 hu) + 2fo 3 (3p 4 frii + 2p 3 /ii 2 ) 

\2u A P A h 13 h 33 + 9u 3 P 4 (h 23 h 13 + /133/112) + 6u 2 P 4 (/i 2 3 + /i 23 /ii 2 + h 33 h n ) 



P^ 



1 4% 

2,u l P 4 {h 13 h 12 + /i 23 /in) + 2h 3 {3P A h 12 + 2P 3 M 



4A§ 



3P 4 /li 3 

"l4 = 



P 



22 



2/13 ' 

6n 4 P 4 /t| 3 + 9u 3 Pih 22 h 23 + 3u 2 P 4 (/t| 2 + 2h 23 h 12 ) + g^AP^ig^gg + 2h 3 {3P 4 h 12 + P 3 /i 22 ) 

2A§ 



12u 4 P 4 /l 23 /l33 + 9u 3 P 4 (/l 2 3 + /l 22 /l33) + 6n 2 P 4 (/l 2 3(/l 22 + /ll3) + /l33/ll2) , 

P23 " 4h£ + 

3u x P A {h X3 h 22 + /i 23 /i 12 ) + 2h 3 (3P 4 h 13 + 2P 3 /i 23 + 3P 4 /i 22 ) 



3P 4 /i 23 

P24 — 



P 
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2/13 ' 

6u 4 P 4 /4 + 3P 4 (3u 3 /i 33 /i23 + ^ 2 (/i| 3 + 2/133/113) + n 1 /i2 3 /H3) + 2/i 3 (3P 4 /i 23 + P3/133) 



3P4/133 
P 44 = 0. 

Upon ensuring the consistency of this system we obtain the set of relations (|29p . We believe 
that the same process applied to any higher order conservation law 

P(u\...,u m ) t = Q(u\...,u m+1 ) x 
will results in the same relations (1291). 
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